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The Hindmash-Rose-model of neuronal 
bursting 


Beginning in the early 20th century, the behaviour of neurons has been described 
by increasingly realistic mathematical models. The very first of these models is called 
integrate-and-fireand is due to LOUS LAPICQUE who developed it in 1907. A better 
model was developed in the early 1960s by RICHARD FITZHUGH and J. NAGUMO 
and is described by 

, v? 
0 ada and 
TW =v+a+—bw. 


This model is still pretty simple and it can be shown that is basically equivalent to the 
VAN DER POL equation 
g+u(y?-Yyt+y=o, 

which was devised by BALTHASAR VAN DER POL in 1920 as a result of his pioneering 
work on vacuum tubes and oscillator circuits. Accordingly, the FITZHUGH-NAGUMO 
model is basically a relaxation oscillator? controlled by an external stimulus Je,¢. 

A much more interesting model is due to HINDMARSH and ROSE? and consists of 
three coupled differential equations 


gS sor fie py =F ls (Ay 
y=—du*+c-y (2) 
zZ=r(s(@ —2,) — z) (3) 





with the parameters a = 1,b = 3, c=1,d=5,r=10°?,s=4,2, = -3 and initial 
conditions of 2 for all three integrators in the final setup. 








Mn contrast to a harmonic oscillator which is typically based on an amplifier with suitable feedback, running in 
resonance mode, a relaxation oscillator switches abruptly between charge and discharge mode and thus yields non- 
harmonic output signals. 

?See [HINDMARSH et al. 1982] and [HINDMARSH et al. 1984]. 
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Figure 1: Numerical simulation of the three coupled differential equations (1), (2), and (3) 


As quick numerical simulation shows, see figure 1, the system must be scaled 
before programmed on an analog computer. Scaling is most easily done manually by 
first applying proper scaling factors to the variables x, y, and z. If, e.g., z is to be 
scaled by a factor of 4, every input of the integrator yielding —z is scaled down by 
that factor. To compensate for this, —z (or z) must be scaled up by a corresponding 
factor of 2 at all inputs of computing elements using this variable etc. In the end, 
these various factors scaling values down and up tend to cancel out in most cases, 
so that typically only a few additional potentiometers are required for a scaled analog 
computer program as compared with an unscaled program. 

The resulting scaled computer setup is shown in figure 2. The scaled parameters 
are these: o* = 4, 6* = 6, c* = 0.066, d* = 1.333; r = 1077, 2 =4, ct = 0.8 with 
initial conditions of +1 accordingly. A bit problematic here are the values r and the 
resulting value of the product rsa* = 0.0032. It is not practical to setup values that 
small directly on a manual potentiometer. In cases like this, the time-scale factor of 
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Figure 2: Scaled analog computer setup for the HINDMARSH-ROSE model 




















the affected integrator can be used to gain an additional scaling factor of 10” as shown 
in the setup. Furthermore, the products rs = 0.004 and rsa* = 0.0032 are setup on 
additional potentiometers. 

Figure 3 shows a typical result obtained by this analog computer circuit with Jeg = 
1. The yellow trace depicts the output potential of the neuron (x), the green trace 
corresponds to the amount of potassium channels (y), and the orange channel shows 
can 
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Figure 3: Typical result of spiking neuron simulation 
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